Trapping Surface Electrons on Graphene Layers and Islands 
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Abstract 

We report the use of time- and angle-resolved two-photon photoemission to map the bound, 
unoccupied electronic structure of the weakly coupled graphene/Ir(lll) system. The energy, dis- 
persion, and lifetime of the lowest three image-potential states are measured. In addition, the 
weak interaction between Ir and graphene permits observation of resonant transitions from an 
unquenched Shockley-type surface state of the Ir substrate to graphene/Ir image-potential states. 
The image-potential-state lifetimes are comparable to those of mid-gap clean metal surfaces. Evi- 
dence of localization of the excited electrons on single-atom-layer graphene islands is provided by 
coverage-dependent measurements. 

PACS numbers: 73.22.Pr, 79.60.Dp, 73.20.-r, 79.20.Ws 
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Graphene on metal surfaces are a materials system of enormous fundamental and applied 
interest. The graphene /metal interface is encountered in the rapidly expanding technological 
system of CVD graphene on Cu foil, in the structurally precise monolayer epitaxial systems 
of graphene on single-crystal Ru, Ir, or Ni, and finally in the metal contacts of graphene 
field-effect transistors or other devices. Questions then arise on the electronic structure 
of graphene on metal surfaces and in fact several recent studies have addressed questions 
such as the role of lattice mismatch on band structures. Most studies of the electronic 
structure of graphene have focused on the band structure in the vicinity of its K point, near 
the Fermi edge. Further there has been a paucity of measurements about its unoccupied 
electronic structure and the dynamics of strongly excited electrons. Image-potential states 
offer one important approach to probe the excited state manifold and are known to vary 
with interfacial quality, dielectric properties, and electronic structure. In graphene the large 
band gap at the T point results in Bragg reflection from the crystal within a certain range of 
energy and momentum. In fact a recent theoretical study has shown the existence of a dual 
Rydberg-like series of even and odd symmetry image-potential states in a single free-standing 
sheet of graphene [l| • Image-potential states on graphene may experience different dynamic 
constraints. For example, the different phase space for decay in two dimensions compared 
to three dimensions may affect the lifetimes for electrons trapped in image-potential states 
on graphene. In fact, more generally the response of the image electron to the composite 
dielectric/metal systems is itself of basic physics interest. 

In this Letter, we investigate the uncharted region of the bound, unoccupied electronic 
structure of epitaxial graphene grown on Ir(lll) in the vicinity of the graphene T point; 
our measurements are made via the image-potential states using angle- and time-resolved 
two-photon photoemission (2PPE) as indicated by arrows in Fig. [U This system was chosen 
for several reasons: First, because of the weak coupling in the graphene/Ir(lll) system, the 
electronic structure of the graphene overlayer is nearly intact, with sharp Dirac dispersion 
characteristics ^]. In addition, the moire corrugation of the epitaxial graphene on Ir(lll) 
has been found to be only 0.35 ±0.10 A based on atomic force microscopy measurements 3] 
indicating a smooth epitaxial graphene surface. Second, the molecular-based growth is well 



characterized and saturates at precisely one monolayer (ML) of epitaxial graphene ^]. Our 
results show that image-potential states may be excited from the Ir /graphene interfacial 
region and have binding energies and lifetimes comparable to those of mid-gap clean metal 




FIG. 1. (Color online). Arrows indicate 2PPE transitions between surface and image-potential 
states. The experimental results (dots) are compared to calculations (lines). The projected bulk- 
band structure of Ir(lll) along the TK direction is shaded according to the total and sp-density 
of states (DOS) at the right and left, respectively. 



surfaces. In addition, spectral measurements of binding energy versus coverage show clearly 
that at low graphene coverage, image-potential electrons are trapped on graphene islands by 
surface work function differences between the metal and graphene regions, an observation 
of high importance for understanding of transport at graphene- metal interfaces 

Our choice of two-photon photoemission is the result of its high temporal and energy 
resolution. Other experimental observations of image-potential states have used scanning 
tunneling spectroscopy, i. e., graphene on SiC jfj| and on Ru(0001) Q. This technique, 
however, measures the image-potential series in the presence of strongly distorting electric 



field between tip and sample and without time-resolved possibilities. 

The experiments were conducted using monochromatic and bichromatic 2PPE, and angle- 
resolved photoemission (ARPES). Details of the monochromatic 2PPE setup at Columbia 
which was used in the photon energy range of 3.8 < hv < 4.9 eV are given in Ref. [8|. Bichro- 
matic and time-resolved 2PPE measurements were performed in Erlangen using pump-probe 
methods with the third harmonic (UV) and the fundamental (IR, 1.51 < hv < 1.62 eV) as 
described in Ref. [9|. Additionally, occupied-state ARPES measurements were performed at 
APE (ELETTRA) using a photon energy of 55 eV with an energy resolution of 20 meV. The 
resolution of the 2PPE experiments was 40 meV. The base pressure in all three ultrahigh 
vacuum (UHV) systems was better than 1 x 10~ 8 Pa. All measurements used p-polarized 
beams. 

Graphene was prepared by cycles of temperature programmed growth, TPG (room tem- 
perature ethene exposure 6 x 10~ 6 Pa for 60 s and flashed to ~ 1450 K), followed by a 



chemical vapor deposition run 



(6 x 10" 6 Pa of ethene for 300 s at 1150 K), to form ex- 



actly one graphene monolayer [4] . Growth was monitored by low-energy electron diffraction 
(LEED) after each cycle, which showed the development of the characteristic moire pattern 



of uniformly oriented graphene (lOj, as graphene coverage varied from to 1 ML; LEED 
patterns (not shown) revealed these patterns clearly. 

Figure |2^i shows the measured 2PPE intensity obtained at 1 ML and for hv = 1.59 eV 
along the TK direction. Three unoccupied bands are observed. The pumping process could 
be deduced from its photon-energy dependence, thus in bichromatic case, all peak positions 
shifted linearly with IR photon energy, indicating that the process involves pumping by a UV 



photon and photoemission by an IR photon All 2PPE features vanished when the IR 



beam was switched to s-polarization, indicating the expected symmetry for image-potential 
states. The effective masses of all three states are 0.9 ±0.1 m e . The binding energies of the 
three states with respect to the vacuum level are given in Table HI The measured energies 
and effective masses are close to the free electron mass and fit well to a Rydberg-like-series 



of image-potential states with a nonvanishing quantum defect 12]. 



Figure [2^l shows that the n = 1 band is most intense for parallel momenta ku between 
points in Fig. Efe] , with the intensity typically decreasing monotonically 



0.08 and 0.17 A" 1 [c 
with increasing k\\ 

resonances [8|. In order to identify possible initial states for 2PPE, we have performed 



131 ] . Direct transitions from initial surface bands can lead to intensity 
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FIG. 2. (Color online), a) Intensity map of the 2PPE signal recorded with photon energy hu = 
1.59 eV for 1 ML graphene on Ir(lll). Points represent the intensity of the lowest n = 1 band, b) 
ARPES map showing initial states for hu = 55 eV. 



ARPES measurements of graphene on Ir(lll). The ARPES data in Fig. [2b show two 
parabolic-like dispersions with a downward curvature. The two branches are shifted from 
fcii = by ±0.033 ± 0.001 A" 1 and have a maximum energy of —0.19 ± 0.01 eV. Similar 
results were also obtained with the fourth harmonic (6.2 eV) in the 2PPE setup. Rashba- 
type splittings of similar magnitude are found in other systems, e.g., a Bi/Ag(lll) surface 
alloy [l4j . These bands are also observed on clean Ir(lll) 15|, [l6| indicating that this surface 
feature is inherent to the clean metal surface. The surface-state energy reported for clean 
and graphene-covered surfaces differs by about 0.2 eV, an effect, which is consistent with a 
charge transfer between substrate and overlayer and which can shift graphene or iridium 
states QQ. 

The initial band dispersion is plotted together with the measured dispersion of the image- 
potential states in Fig. [T] (blue dots). The arrows connecting initial states to the n = 1 
image-potential band are at slightly larger kn values than the enhanced intensity in Fig. [21a. 
In the absence of resonances the 2PPE intensity along image-potential bands decreases 
continuously with increasing parallel momentum 13j. In the present case due to the finite 
energy and angle resolution the intensity maxima are shifted to lower kn values compared 



TABLE I. Experimental and calculated binding energies and lifetimes for image-potential states 
on graphene/Ir(lll). 



n E^v (eV) E^ alc (eV) r (fs) 

1 0.83 ±0.02 0.59 35 ± 3 

2 0.19 ±0.02 0.18 114 ±6 

3 0.09 ±0.02 0.08 270 ± 12 



to the position found in the dispersion analysis. The additional resonance into the n = 2 
band (see Fig. [1]) can be inferred from the similar intensity as for the n — 1 state at k\\ = 
in Fig. [2^ and is confirmed by photon-energy-dependent data presented in the supplemental 



material 19] 



In order to understand the character of the initial state, we calculated the projected bulk- 



band structure of Ir(lll) using a non-relativistic parameterized tight-binding scheme [20 ]. 
Figure [1] shows this projected structure along the TK direction, at the right. The shading 
represents the one- dimensional density of states (ID-DOS). The left part of Fig. [Ushows the 
ID-DOS of bands according to their sp-character. For kn = 0, the lower edge of the sp-band 
gap is at —0.7 eV, which corresponds to the L 2 ' point. The band edge of the total projected 
bulk-band structure disperses upward from the L 3 point around ±0.8 eV and picks up sp- 
contributions. On the other hand the lower sp-band edge shows a downward dispersion. 
The energy of the Shockley surface state was calculated using the sp-band edges within a 
scattering model 2l|. The calculations used the experimental work function of 4.65 eV for 
graphene on Ir(lll). The calculated bands were shifted by ±0.033 A" 1 to account for the 
experimentally observed Rashba splitting and are drawn as green lines in Fig. [1] in the region 
below the Fermi energy. The experimentally extracted dispersion shown by dots agrees well 
with the calculation. The Shockley-type surface state is apparently not quenched by the 



graphene layer at a distance of 3.4 A 22], because its probability density is concentrated at 
the Ir(lll) surface. 

The scattering model was also used to calculate the energies of the image-potential- 



band series [21]. The calculated binding energies, given in Table HI are approximately those 
expected for states located near midgap (see Fig. [[]). However, the calculated n = 1 binding 
energy is significantly smaller than the experimental value. This discrepancy is due to the 
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FIG. 3. a) Time-resolved measurements of the image-potential-state series at k\\ = 0. b) Lifetimes 
as function of binding energy compared to various power-law dependencies (see text). 



fact that the scattering model calculation neglects the round-trip phase shift 2<p gr of the 
graphene layer. Using the expressions for the phase shift at the substrate and the image- 



potential barrier 23[ we obtain (ftc — 0.637T and (pB = 1.027T, respectively. The total phase 
shift for the n = 1 state is 2n, from which we obtain (j) gr = 0.187T. Note that such a small 
phase shift leads to a significant change in binding energy from 0.59 eV to 0.83 eV (see 
Table!]). 

The time-resolved spectra of the image-potential states were also measured and are shown 
in Fig. As summarized in Table HJ lifetimes of tens to hundreds of femtoseconds are 
obtained. These are comparable to values obtained for Cu(100) with a similar midgap 
image-potential state position and hence bulk evanescent decay length in the metal crystal 



Note, as an aside, that the curve measured at the energy of the n = 4 image-potential 
state in Fig. [3^l shows weak quantum beats [ill] for delay times < 300 fs. The data in Table U 
show that lifetimes vary with binding energy approximately oc E~ l (solid line in Fig. [3b). 

rr 

The asymptotic, classical r oc E~ 3 / 2 behavior [24j (dashed line in Fig. [3b) is not reached for 



n < 4. Similar behavior has been found on copper surfaces 



25]. 



An important issue for carrier movement at graphene/metal interfaces is the degree of 
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FIG. 4. Sample work function (open symbols) and image-potential states n = 1, 2, and 3 binding 
energies (filled symbols) as a function of graphene coverage. Dashed line represents a linear fit for 
the work-function change. 



lateral confinement. This confinement can be examined at low graphene coverage, obtained 
via a small number of sequential TPG cycles. From previous studies, it is known that one 
TPG cycle covers a fraction of about 20% of uncovered Ir surface [17j . After one TPG eye 



the typical island size is (35 nm) 2 and after the second cycle of the order of (100 nm) 2 17]. 
As has been shown in earlier work 



26], the average and local work function play an 
important role in interfacial electron localization. Thus the average work function $ was 
measured via monochromatic 2PPE and the expression $ = 2hu — AE, i.e., where hu is the 
photon energy and AE is the difference between the Fermi level cutoff and the low-energy 
cutoff. Figure H] displays the work function (open symbols) as a function of graphene cov- 
erage. The work function decreases approximately linearly from a value 5.79 ± 0.10 eV to 
4.65 ± 0.10 eV from Ir(lll) to 1 ML graphene. Reported values of the work function for 



Ir(lll) are 5.76 and 5.79 eV [27] . The work function of the graphene-covered surface on 
r(lll) is between the values for Pt(lll) of 4.87 eV and free-standing graphene of 4.48 eV 
281 ] . which is consistent with the weak bonding between the Ir(lll) and the graphene over- 
layer and a p-doping of the graphene [2| . The linear decrease of the work function is known 



for other systems and is due to the averaging over substrate and overlayer islands 



26] 



Image-potential states were observed at all coverages reported here using 2PPE. However, 
for the clean surface or uncovered substrate areas the available photon energies were not 
sufficient to populate image-potential states due to the large work function of Ir(lll). The 
image-potential-state energies, measured relative to the Fermi level, are shown in Fig. HI 
The energies are generally constant over the coverage range from 0.2 to 1 ML with the 
intensity increasing monotonically with coverage. Note that the graphene Dirac cone at 
the K point has been clearly observed for more than three TPG cycles or 0.5 ML graphene 
coverage [l^]. The constant energy of the image-potential series as a function of coverage in 



Fig. H]is a direct result of the localization of the electrons on the graphene islands 26|]. The 
electrons respond to the local work function if the average island dimensions are larger than 
the typical distance of the probability density maximum which is of the order of nanometers 
for the lowest n image-potential states. Note that the localization on the graphene islands is 
facilitated by the large work function difference between the graphene layers and the Ir(lll) 
substrate. For small graphene islands, an energy shift proportional to d~ 2 , where d denotes 
the characteristic island size, is exp ected due to the lateral localization of the electron in a 



two-dimensional quantum well 291. However, these shifts would be < 1 meV for the island 

n 

sizes expected for the current preparation conditions |17j. 

In summary, we have observed and measured the properties of image-potential states 
on a graphene monolayer on Ir(lll). The binding energy of the n = 1 image-potential 
state is 40% larger than expected from the position of the graphene vacuum level rela- 
tive to the Ir(lll) band gap. There is no prominent indication of a second main series of 
image-potential states as predicted for free-standing graphene jlj]. Apparently, the under- 
lying metal substrate breaks the mirror symmetry of the graphene layer and the state of 
odd symmetry shifts up in energy as has been calculated for graphene on Ru(0001) [7]. In 
addition, the image-potential states can be excited efficiently from a downward dispersing 
Shockley surface state in the sp-band gap of the Ir(lll) band structure indicating a sizable 
overlap of the wave functions of these states located at the substrate interface and graphene 
surface, respectively. The measured lifetimes of the image-potential states are comparable to 
similar clean metal surfaces. Recently, similar results have also been obtained for graphene 



on Pt(lll) 30]. Apparently, the evanescent coupling of the image-potential-state wavefunc- 
tions to the underlying electronic states of the Ir(lll) bulk and surface states is not altered 
by the graphene layer. Three-dimensional localization of electrons on graphene islands has 



been observed for submonolayer coverages obtained by individual TPG cycles. However, 
even for the smallest island size, no energy shift due to localization was observed within the 
experimental uncertainty. Further development is needed to prepare well-ordered graphene 
islands with controlled lateral extension. A different approach would be to exploit the moire 
pattern on more corrugated graphene layers 
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